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Abstract. An algorithmic proof of General Neron Desingularization is given 
here for one dimensional local domains and it is implemented in Singular. Also 
a theorem recalling Greenberg’ strong approximation theorem is presented for one 
dimensional Cohen-Macaulay local rings. 
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Introduction 

A ring morphism u : A ^ A' has regular fibers if for all prime ideals P G Spec A 
the ring A'/PA' is a regular ring, i.e. its localizations are regular local rings. It 
has geometrically regular fibers if for all prime ideals P G Spec A and all hnite held 
extensions K of the fraction held of AjP the ring K ®a/p A'/PA' is regular. If for 
all P G Spec A the fraction held of A/P has characteristic 0 then the regular hbers 
of u are geometrically regular hbers. A hat morphism u is regular if its hbers are 
geometrically regular. 

In Artin approximation theory [2] an important result is the following theorem 
generalizing the Neron Desingularization 0. H- 

Theorem 1 (General Neron Desingularization, Popescu [8], |9], [10], Andre |T], 
Swan |T3|, Spivakovski IE])- Let u \ A ^ A' be a regular morphism of Noetherian 
rings and B a finite type A-algebra. Then any A-morphism v \ B ^ A' factors 
through a smooth A-algebra C, that is v is a composite A-morphism B ^ C ^ A'. 

The purpose of this paper is to give an algorithmic proof of the above theorem 
when A, A' are one dimensional local domains and A D Q. This proof is somehow 
presented by the second author in a lecture given in the frame of a semester of 
Artin Approximation and Singularity Theory organized in 2015 by CIRM in Lu- 
miny (see http://hlombardi.free.fr/Popescu-Luminy2015.pdf). The algorithm was 
implemented by the authors in the Computer Algebra system Singular |3] and 
will be as soon as possible found in a development version as the library GND.lib at 
https://github.com/Singular/Source. 

The support from the Department of Mathematics of the University of Kaiserslautern of the first 
author and the support from the project ID-PGE-2011-1023, granted by the Romanian National 
Authority for Scientific Research, GNGS - UEEISCDI of the second author are gratefully acknowl¬ 
edged. Both authors thank GIRM, Luminy who provided excellent conditions and stimulative 
atmosphere in the main stage of our work. 
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We may take the same General Neron Desingularization for v,v' : B ^ A' if 
they are closed enough as Examples IH [10] show. The last section computes the 
General Neron Desingularization in several examples. We should point that the 
General Neron Desingularization is not unique and it is better to speak above about 
a General Neron Desingularization. 

When A is the completion of a Gohen-Macaulay local ring A of dimension one 
we show that we may have a linear Artin function as it happens in the Greenberg’s 
case (see 0). More precisely, the Artin function is given by c —)■ 2e + c, where e 
depends from the polynomial system of equations dehning B (see Theorem l20l) . 

1. The theorem 

Let u : A A' he a flat morphism of Noetherian local domains of dimension 1. 
Suppose that A D Q and the maximal ideal m of A generates the maximal ideal of 
A. Then u is regular morphism. Moreover, we suppose that there exist canonical 
inclusions k = A/m A, k' = A/mA —>■ A such that u{k) C k'. 

Let B = A\Y]/I, Y = (W,..., W)- If / = (/i, • • •, /r), r < n is a system of 
polynomials from I then we can dehne the ideal Aj generated by all r x r-minors 

of the Jacobian matrix • After Elkik |1] let Hb/a be the radical of the ideal 

((/) : l)AfB, where the sum is taken over all systems of polynomials / from 
I with r < n. Then Bp, P G SpecU is essentially smooth over A if and only if 
P Hp/A by the Jacobian criterion for smoothness. Thus Hp/A measures the non 
smooth locus of B over A. 

Definition 2. B is standard smooth over A if there exists / in J as above such 
that 1 e ((/) : l)AfB. 

The aim of this paper is to give an easy algorithmic proof of the following theorem. 

Theorem 3. Any A-morphism v : B ^ A factors through a standard smooth 
A-algebra B'. 

If A is essentially of hnite type over Q, then the ideal Hp/A can be computed 
in Singular by following its dehnition but it is easier to describe only the ideal 
((/) • dehned above. This is the case considered in our algorithmic part, 

let us say A = k[x]/F for some variables x = {xi,.. .Xt), and the completion of A 
is K\x\/{F) for some held extension k <Z K. When v is dehned by polynomials y 
from K[x\ then our problem is easy. Let L be the held obtained by adjoining to k 
all coefficients of y. Then R = L[x]/{F) is a subring of A containing Imn which is 
essentially smooth over A. Then we may take B' as a standard smooth A-algebra 
such that i? is a localization of B'. Gonsequently we suppose usually that y is not 
polynomial dehned and moreover L is not a hnite type held extension of k. 

2. Reduction to the case when LTb/a n A 7^ 0. 

We may suppose that v{Hb/a) 7^ 0. Indeed, if v{Hp/a) = 0 then v induces an A- 
morphism v' : B' = B/Hp^A A and we may change {B,v) by {B',v'). Applying 
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this trick several times we reduce to the case v{Hb/a) 7^ 0 . However the fraction 
field of Imu is essentially smooth over A by separability, that is Hi^-^/aA' 7^ 0 and 
in the worst case our trick will change H by Imn after several steps. 

Choose P' G (Aj((/) : /)) \ / for some system of polynomials / = (/i ,..., fr) 
from / and d' G (n(P')H') O A, d' 7^ 0 . Moreover we may choose P' to be from 

M[{f) : l) where M is a r x r-minor of • Then d' = v{P')z G [v{Hb/a)) A A 

for some x G H'. Set Bi = B[Z]/{fr+i), where fr+i = —d'+P'Z and let ui : Pi —)■ A' 
be the map of P-algebras given hj Z ^ z. It follows that d' G ((/, fr+i) ■ {I, fr+i)) 
and d' G A/, d' G A^fr+i- Then d = d'^ = P modulo (/, fr+i) for P = P'‘^Z^ G PIbi/a- 
For the reduction change B by Pi and the Jacobian matrix J = [df/dY) will be 

^ ^ Note that d G Hb/a H A. 


now the new J given by 


P' 


Example 4. Let 01,02 G C be two elements algebraically independent over Q and 
p a root of the polynomial + T + 1 in C. Then k' = = Q(p, oi, 02). 

(o| + O3 + 1) 


Let.4= 

V(a^I -2^2) 

defined as 


1 and P = 


{xi,X2) 

(U - U) ’ 

V : B - 

- 

— 

-s- 01X2 

- 

-^ 0103X2 


30 j 

F3I-> 

/ vj—h CI2X2 


k'lxi,X2j 

and the map v 




Xn 


50 


i=0 


i=31 


^ i\ 


This is an easy example. Indeed, let v" : B” = H[o3, 01X2, ^(13)] —)■ A' be the 


inclusion. We have Imn C B" = 


A[T,Yj,Y 3 


and P^; 3+1 ^ 


A[T,Y^,Ys] 
(T2 + T + 1) 


2T+1 


(T2 + T 

is a smooth H-algebra, which could be taken as a General Neron Desingulariza- 
tion of P. Applying our algorithm we will get more complicated General Neron 
Desingularizations but useful for an illustration of our construction. 

Then Imn, the new P will be —-, where the kernel is generated by the following 

Kern 

polynomial: 

ker[1]=Y1~2+Y1*Y2+Y2"2 

Next we choose / = y|^+Yil2+T2^ and we have M = 2l2+hi and 1 G ((/) : /) and 
hence P' = Yi + 2y2- Therefore n(P') = (20103 + 01) • 0:2 and d' = X2, z = 


Therefore d = d'"^ = 


2O1O3 + Oi 


XX 


To be able to construct 


, [xl in Singular we will add a new variable 

2O1O3 + Oi 

o and we will factorize with the corresponding polynomial 2oi03-o+Oi -o — 1 . We will 
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B[Yt] 


see this a as a new parameter from k' G A'. Then we change B hy = , 

{-d' + P'Ki) 

and extend n to a map ni : i?i —)■ A' given by 14 —)■ a. Changing B by Bi we may 
assume that d E Hb/aA A. 

Example 5. Note that we could use B instead Imn. In this case we choose / = 
— Y2 and take M = SY^ and 1 G ((/) ; /). Therefore we obtain P' = SY^, 
d' = d = xt and the next computations are harder as we will see in the Examples 

[IT] and 123 

C|a:i,X2] . 


Remark 6. We would like to work above with A" = 


instead A\ v being 


- ^ 2 , 

given by ^( 12 ) = o>ipx2- But this is hard since we cannot work in Singular 
with an inhnite set of parameters. We have two choices. If the dehnition of v 
involves only a hnite set of parameters then we proceed as Example HI using some 
A' D Imn. Otherwise, we will see later that in the computation of the General Neron 
Desingularization we may use only a hnite number of the coefficients of the formal 
power series dehning v{Y) and so this computation works in Singular. 

Remark 7 . As we may see our algorithm could go also when A' is not a domain, 
but there exist P G M{{f) : I) as above and a regular element d E m with d = 
P modulo /. If A is Cohen-Macaulay we may reduce to the case when there exists 
a regular element d E Hb/aA A. However, it is hard usually to reduce to the case 
when d = P modulo I for some P E M{{f) : I). Sometimes this is possible as shows 
the following example. 

Example 8. Let oi, 02 G C be two elements algebraically independent over Q. Con- 


. I Q[a;i,a;2,a;3] 

sider A = ^^, 

\[X2 / (xi,X2,X3) 

K'lxi,X2,X3j 


and B = 


(113 _ y|) ’ 


K' = 


Q(ai,a 2 )[a 3 ] 

(a| - 0102) 


A' = 


Xn 


X 3 ,X-( 


4 ) 


and the map v dehned as 


V : 


B 


A' 


I2L 


03X1 

03X2 


Y 3 h 


30 


50 


2=0 2=31 


Then Im v, the new B, will be —-, where the kernel is generated by six 

Ker V 

polynomials: 

ker[l]=x2*Yl-xl*Y2 

ker [2] =Y1- 3 ^ 2-3 

ker [3] =xl*Yl''2-x2*Y2''2 

ker [4] =xl''2*Yl-x2''2*Y2 

ker [5] =xl*x2''2*Y2-x3''2*Yl 

ker [6] =xl''2*x2*Y2''2-x3''2*Yl''2 
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Next we choose / = 0:2^1 — a;il2 and we have M = —xi. We may take N = — G 
((/) : /) and P' = xixl- Note that xi — 0:2 is a zero divisor in A bnt d' = P' is 
regnlar in A. In this example we may take d = d' = P' = P. 

Remark 9. Changing B by Imn can be a hard goal if let us say A' is a factor 
of the power series ring over C in some variables x and v{Y) is dehned by formal 
power series whose coefficients form an inhnite held extension F of Q. If v{Y) are 
polynomials in x as in Examples 01 |H] then it is trivial to hnd the General Neron 
Desingularization of B as we explained already in the last sentences of Section 1 . For 
instance in Example |8], B' could be a localization of K' ®q A. Thus Examples IU [8] 
have no real importance, they being useful only for an illustration of our algorithm. 
This is the reason that in the next examples the held L obtained by adjoining to k 
of all coefficients of y will be an inhnite type held extension of k and v{Y) are not 
all polynomials in x. 

However, this will complicate the algorithm because we are not able to tell to the 
computer who is v{Y) and so how to get d'. We may choose an element a E m 
and hnd a minimal c G M such that G (n(M)) + (this is possible because 
diniH = 1 ). Set d' = Y. It follows that d' G (n(M)) + (d'^) C {v{M)) + (d'^) C ... 
and so d' G (n(M)), that is d' = v{M)z for some 2; G W. Certainly we cannot hnd 
precisely z but later it is enough to know just a kind of truncation of it modulo d'®. 


Example 10. Let a* G C, z G M be elements algebraically independent over 


and p a root of the polynomial + T + 1 in C. Let A = 


B = 


{Y^ + YAY2 + Yi) 


, A' = cmd the map v dehned as 


Q[a;i,a;2] 

Yl ‘^ 2 ) / {xi,X2) 


and 




Xn 


V : 


B -- W 



2=0 2=10 


As in Example m we may take d' = X2, d = d'^ and a. Our algorithm goes exactly 
as in Examples 01 [T6l 1211 providing the same General Neron Desingularization. This 
time we cannot hnd an easy General Neron Desingularization as in the hrst part of 
Example 01 
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Example 11 . Let A = Then the inelueion A c A' 

[xi - x'i) [xi - xi) 

oo oo 

is regular. Let di = aijX^ + xi Ai^2 ^ C|a:i, 0:2] for i = 3,4 with aiQ = 1 and 
i=o i=o 

yi = 1/2 = 2/3 = X29z, 2/4 = a;26'4. Let /i = - a;^yiF2, /2 = Y^ - X2Y2Y3 be 

polynomials in A[Y], Y = (Fi,..., F4) and set B = A[Y]/{f), f = (/i,/2). 

If i? is a domain and m G i? is such that Y"^ — u G R[Y] has no solutions in 
Q{R) then it is easy to see that R[Y]/ (F^ — u) is a domain too. In our case we get 
that R = 74[Fi, F2, F3]/(/i) and B = i?[F4]/(/2) are domains too. Then the map 
V : B ^ A' given by F —)■?/ = (j/i,..., 2/4) is injective if we suppose that 6*3, 9 ^ are 
algebraically independent over A. This follows since i? is a domain and dimi? = 
tr degQ(^) Q{B) = tr degg^^^ Q(Imu) = 2 = dimimu. Moreover we will assume that 
the fields Li = Q((ap-, Ai)i) F = 3,4 have inhnite transcendental degree over Q. The 

Jacobian matrix ( ) have a 2 x 2 —minor M = det ( Tr~] = 4F3F4 ^ (/). 

V oY I V oYj I i<i<2 

3 <i <4 

Note that v{M) = X2I/5, where 7/5 = 1/(46*36*4). Then we may take Bi = i?[F5]/(/3), 
/s = —X2 + MF5 and ui given by F5 — >■ y^. Clearly, P = M‘^Y^ G Hb^/a and 
0 7^ d = a;2 = Vi{P) G A. 

3 . Proof of the case when Hb/aA A 0 . 

Thus we may suppose that there exists / = (/i,...,/^), r < n a system of 
polynomials from /, a r x r-minor M of the Jacobian matrix {dfi/dYj) and N G 
((/) : I) such that 0 7^ d = MN modulo /. Set A = A/(d^), A' = A'/(RA'^ 
u = A (8)^ u, B = B/d^B, v = A 0^1 v. Clearly, u is a regular morphism of Artinian 
local rings. 

Remark 12. The whole proof could work with A = A/d?u for any u & m. We 
prefer to take u = d as is done in [ 9 ] and m but we could choose u d, u E m \ m? 
for easy computations. 

By P 19 , 7 . 1 . 5 ] for every held extension L/k there exists a hat complete Noetherian 
local A-algebra A, unique up to an isomorphism, such that mA is the maximal ideal 
of A and AjmA = L. It follows that A is Artinian. On the other hand, we may 
consider the localization A^ of L A in m{L A) which is Artinian and so 
complete. By uniqueness we see that Ab = A. Set k' = A'/mA'. It follows that 
A = Ay. Note that A^ is essentially smooth over A by base change and A' is a 
hltered union of sub-A-algebras Al with L/k hnite type held sub extensions of k'/k. 

Let V be given by F -A- y E A'”. Choose L/k a, hnite type held extension such 
that Al contains the residue class y E A'" induced by y. In fact y is a vector of 
polynomials in the generators of m with the coefficients in k’ and we may take 
L = k{{cy)y). Then v factors through Al- Assume that k[{ci,)i/\ = k[{Uy)i/\/J for 
some new variables Lf and a prime ideal J C k\U]. We have HL/k 7^ 0 because L/k 
is separable. Then we may assume that there exist u = (cui,..., Up) in J such that 
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p = det{duji/dUi,)i^y^[p\ ^ 0 and a nonzero polynomial r G ((ca) : J)\ J. Thus L 
is a fraction ring of the smooth fc-algebra {k[U]/{u;))pT-. Note that u:,p,T can be 
considered in A because k C A and Cj, G A! because k' G A'. 

Then v factors through a smooth A-algebra C = {A\U]/{u)) for some polyno¬ 
mial 7 which is not in m{A\U]/{u:))pr. 

Lemma 13 . There exists a smooth A-algebra D such that v factors through D = 

D- 

Proof. By our assumptions u{k) C k'. Set D = {A\U]/{uj))pr^ and w ■. D ^ A' he 
the map given by Uy —?• Cy. We have C = A D- Certainly, v factors through 

w = A w but in general v does not factor through w. □ 

Remark 14 . \i A' = A then A = A' and we may take D = A. 

Remark 15 . Suppose that k G A but L A' and so k' A A!. Then D = 

{A[U, Z]/{u — (PZ))pT-.y, Z = [Zy] is a smooth ^l-algebra and D = C[Z]. Since 

V factors through a map C ^ A' given let us say hj U ^ X + df’A' for some A 
in A' we see that a;(A) = 0 modulo d^, that is a;(A) = df'z for some 2; in A'. Let 
w \ D ^ A' he the y 4 -morphism given by {U,Z) —)■ (A, 2;). Certainly, v factors 
through w = AZ)aw but in general v does not factor through w. If also k gl A then 
the construction of D goes as above but using a lifting of u, r, 7 from k[U] to A[U]. 
In both cases we may use D as it follows. 

Example 16 . We reconsider Example 01 We already know that d = x\. The 
algorithm gives us the following output: 

This is C: 

// characteristic : 0 

// number of vars : 5 

// block 1 : ordering dp 

// : names al a3 a xl x2 

// block 2 : ordering C 

// quotient ring from ideal 

_[1]=3*al*a+2*a3+l 

_[2]=a3~2+a3+l 

_[3]=xl~3-x2''2 

_[4] =x2~6 

This is D: 


// 

characteristic 

0 


// 

number of vars 

5 


// 

block 1 

ordering 

dp 

// 


names 

al a3 a xl x2 

// 

block 2 

ordering 

C 

// 

quotient ring from ideal 


_[1]=3*al*a+2*a3+l 




[2] =a3~2+a3+l 

[3] =xl~3-x2''2 
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Indeed, 


C = 


A[ai^ Os, a] 


and 


D = 


(3fliO -|- 2(33 “1“ 1) “1“ ®3 1) ^ 2 ) 

74[ai, 03, a] 


(3(3i( 3 + 2a3 + 1, + 03 + 1) 

Note that the hrst polynomial from C comes from the standard basis computation 
of the ideal (203030 + oio — 1, o| + 03 + 1). 


Example 17. Now we reconsider Example O We know that d = x\. The algorithm 
gives us the following output: 


This is C: 



// 

characteristic 

0 

// 

number of vars 

5 

// 

block 

1 

ordering dp 

// 



names al 

// 

block 

2 

ordering C 

// 

quotient ring 

from ideal 


_[l]=a3~2+a3+l 
_[2]=xl~3-x2''2 
_[3]=al~2*a-a3 
_[4]=x2~12 


This is D: 


// 

characteristic 

0 


// 

number of vars 

5 


// 

block 

1 

ordering dp 

// 



names 

al 

// 

block 

2 

ordering 

C 

// 

quotient ring 

from ideal 



_[l]=a3~2+a3+l 
_[2]=xl~3-x2''2 
_[3]=al~2*a-a3 


Indeed, 


and 


C 


A[oi, O3, o] 

(o| + 03 + 1, a\a — 03, X2^) 


D 


A[ai^ 03, o] 

(o| + 03 + 1, a\a — 03) 


Example 18. In the case of Example Owe obtain the following output: 
This is C: 

// characteristic : 0 

// number of vars : 5 

// block 1 : ordering dp 









names 


// : 

// block 2 : 

// quotient ring from 
_[l]=x2~3-x3''2 


ordering 

ideal 


_[2]=xl~3-x3''2 
_[3]=x3~8 
This is D: 

// characteristic : 0 

// number of vars : 5 

// block 1 : ordering 

// : names 

// block 2 : ordering 

// quotient ring from ideal 
_[l]=x2~3-x3''2 


[2]=xl~3-x3''2 


al a3 xl x2 x3 
C 


dp 

al a3 xl x2 x3 
C 


Indeed this is the case since we have d = x\xl and hence 


and 


^[01,03] 


D = ^[01,03]. 


Example 19. In Example ITT] we consider 01,02 algebraically independent over Q 

and set 6^3 = 1 + 01X2 and 6^4 = 1 + 02X2. Suppose that 0 - = 6i modulo x^. We have 

^3 q2 ^ 

Vs = X2O3, Va = X2OA, yi = ^, 1/2 = ?/5 = ■ Choose y', i e [ 5 ] polynomials 

O4 O3 

with degrees < 11 in X2 and linear in xi such that ?/' = mod(xi, We get 
y[ = Hi = Q\/Q\ = 6*'3/6*'4mod(xi,X2^) and similarly for ?/', i > 1 . Here we use the 

e 

fact that 6*4^^ = ^^(1 — Oly for some e >> 0 because 1 — 6*| is nilpotent in the ring 
j=i 

H[oi, 02, 03, 04]. Thus the coefficients of j/', i G [ 5 ] belong to the held L obtained by 
adjoining to Q the coefficients of 6*3, 6*4. Note that in this case L = Q (Q[oi, ...,04]). 
Then we obtain the following output: 

This is C: 

// characteristic : 0 

// number of vars : 4 

// block 1 : ordering dp 

// : names al a2 xl x2 

// block 2 : ordering C 

// quotient ring from ideal 

_[l]=x2~3-xl''2 

_[2]=xl~8 

This is D: 

// characteristic : 0 
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// 

number of vars 

4 

// 

block 

1 

ordering dp 

// 



names al 

// 

block 

2 

ordering C 


// quotient ring from ideal 
_[l]=x2~3-xl''2 


Thus c = 

( 4 ) 

with il[ai,..., 04]. 


= A[ai,..., 04] which is smooth over A. Then D is equal 


Back to our proof note that the composite map -B —)■ C* —)■ Z) is given by y —)■ 
y' + d^D for some y' G -D”. Thus I{y') = 0 modulo d^D. Since v factors through 
w we see that w{y' + d^D) = y. Set y = w{y'). We get y — y E d^A'"’, let us say 
y — y = d?e for e e dA''^. 

We have d = P modulo I and so P{y’) = d modulo d^ in D because I{y') = 0 
modulo d^D. Thus P{y') = ds for a certain s E D with s = 1 modulo d. Assume 
that P = NM for some N E ((/) : I). Recall from beginning of Section 2 that 
the new M is now the old one multiplied with P' and the new N is the old one 
multiplied with Z^. Let H be the n x n-matrix obtained adding down to {df/dY) 
as a border the block ( 0 |ld„_r). Let G' be the adjoint matrix of H and G = NG'. 
We have 

GH = HG = NMldn = PId„ 

and so 

dsldn = P{y')Un = G{y')H{y'). 

Then t := H{y')e E dA'"" satishes 

G{y')t = P{y')e = dse 

and so 

-y) = dw{G{y'))t. 

Let 

( 1 ) h = s{Y -y')-dG{y')T, 

where T = (Ti,..., T„) are new variables. The kernel of the map ip \ D\Y,T] ^ A' 
given by y —)■ y, T —)■ f contains h. Since 

s(y — y') = dG{y')T modulo h 

and 

f(Y) - /(»') = 5 ;; §^{y')(y, - v',) 

j ^ 

modulo higher order terms in Yj — y'j by Taylor’s formula we see that for p = 
max* deg fi we have 

(2) s’’f(Y) - 5 5 ;^ + d^Q = 8<‘-'dP{p')T + rfQ 

j ^ 
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modulo h where Q G T^D[T]'’. This is because {df/dY)G = (Pldr| 0 ). We have 
f{y') = (Pb for some b G dD^. Then 

( 3 ) Qi = s^bi + s^Ti + Qi, i e[r] 

is in the kernel of ip because d?ip{g) = d?g{t) G {h{y,t), f{y)) = ( 0 ). Set E = 
D[Y,T]/{I, g, h) and let -0 : P —)■ W be the map induced by ip. Clearly, v factors 

•tp 

through 0 because v is the composed map B ^ B 0^^ D = D[Y ]/1 E ^ A'. 

Note that the r x r-minor s' of {dg / dT) given by the hrst r-variables T is from s'’^+ 
(T) C 1 + (d,T)P)[y',T] because Q G {Tf. Then U = (P)[y,T]/(h,is smooth 
over D. We claim that / C {h^ g)D\Y,T]ss's" for some other s" G 1 + {d^T)D\Y,T]. 
Indeed, we have PI C {h,g)D[Y,T]s and so P{y' + s-HG{y')T)I C (h, c/)P[F, T],. 
Since P{y' + s~^dG{y')T) G P{y') + d{T) we get P{y' + s~^dG{y')T) = ds" for some 
s" G 1 + (d,T)P[F,T], It follows that s"I C {h^g)D\Y^T]ss' because d is regular in 
U, the map D ^ U being flat, and so / C {h, g)D[Y,T]ss's"- Thus Egs's” — Us" is a 
P-algebra which is also standard smooth over D and A. 

As w{s) = 1 modulo d and w{s'),w{s") = 1 modulo (d, f), d,t ^ rnA' we see that 
w{s),w{s'),w{s") are invertible because A' is local and 0 (thus v) factors through 
the standard smooth A-algebra Ess's"- 

4 . A THEOREM OF GREENBERG’S TYPE 

Let (A, m) be a Cohen-Macaulay local ring (for example a reduced ring) of di¬ 
mension one. A' = A the completion of A, P = A\Y]/I, Y = (Yi,..., Y„) a hnite 
type A-algebra and c, e G N. Suppose that there exist / = (/i,...,/^) in J, a 
r X r-minor M of the Jacobian matrix (df/dY), N G ((/) : I) and an A-morphism 
n : P —)■ such that (n(MiV)) D rrf 

Theorem 20 . Then there exists an A-morphism v' \ B ^ A such that v' = 
V modulo m'^, that is v'iY J- /) = v{Y + I) modulo m'^. 

Proof. We note that the proof of Theorem [ 3 ] can work somehow in this case. Let 
y' G A” be an element inducing v{Y + /). Then C {{MN){y')) + C 

((MA^)(?/')) C ... by hypothesis. It follows that C {{MN){y')). Since A 

is Cohen-Macaulay we see that contains a regular element of A and so {MN){y') 
must be regular too. 

Set d = {MN){y'). Next we follow the proof of Theorem E] with P = A, s = 1 , 
P = MN and P, G such that . 

dld„ = P{y')ldn = G{y')H{y'). 

Let 

h = Y-y'-dG{y')T, 

where T = (Ti,..., T„) are new variables. We have 

f(Y)-f{y')^dP{y')T + d^Q 

modulo h where Q G T^A[T]^. But f{y') G C d'^m^A'^ and we get f{y') = 

d% for some b G m'^A''. Set gi = bi + Ti Y Qi, i G [r] and E = A[Y, T]/(/, h, ^f). 

11 


We have an A-morphism f 3 : E ^ Ajrrf given by {Y,T) —)■ (y', 0 ) becanse I{y') = 
0 modnlo h{y', 0) = 0 and 5f(0) = 6 G rrfA^. 

As in the proof of Theorem [ 3 ] we have Esign = Us"-, where U = {A\Y,T]/{g, h))si- 
This isomorphism follows becanse d is regnlar in A and so in U. Conseqnently, 
is smooth over A. Note that /3 extends to a map f 3 ' : Eg's" —t Ajrrf. By the Implicit 
Fnnction Theorem jd' can be lifted to a map w : Eg's" A which coincides with ( 3 ' 
modnlo It follows that the composite map n', B —)• Es's" ^ A works. □ 

Corollary 21. In the assumptions of the above theorem, suppose that {A,m) is 
exeellent Henselian. Then there exists an A-morphism v" : B ^ A sueh that v" = 
V modulo m'^, that is v"{Y + J) = v{Y + I) modulo m'^. 

Proof. An excellent Henselian local ring {A, m) has the property of Artin approx¬ 
imation by [ 9 ], that is the solntions in A of a system of polynomial eqnations / 
over A are dense in the set of the solntions of / in A. By Theorem [ 20 ] we get an 
A-morphism v' : B ^ A snch that v' = v modnlo Then there exists an A- 

morphism v" \ B ^ A snch that v" = v' = v modnlo by the property of Artin 
approximation. □ 

Theorem 22. Let (A, m) be a Cohen-Macaulay local ring of dimension one, B = 
A\Y]/I, Y = {Yi,... ,Yn) a finite type A-algebra, e G N and f = (/i,...,/^) a 
system of polynomials from I. Suppose that A is excellent Henselian and there exist 
a r X r-minor M of the Jacobian matrix {df/dY), N G ((/) : I) and y' G A” such 
that I{y') = 0 modulo and {{NM){y')) D m^. Then the following statements are 
equivalent: 

( 1 ) there exists y" G A" such that I{y") = 0 modulo and y” = y' modulo m®, 

(2) there exists y ^ such that I{y) = 0 and y = y' modulo m®. 

For the proof apply the above corollary and Theorem [201 


5. Computation of the General Neron Desingularization in 

Examples lillSlISlfTT] 

Example 23. We would like to compute Example H] in Singular using GND.lib. 
We quickly recall the example. 

Let 01,02 G C be two elements algebraically independent over Q and p a root 


of the polynomial -|- T -|- 1 in C. Then k' = 


/ Q[a;i,a;2] 
dehned as 


Q(0i,02)[03] 

(o| -|- O3 -|- 1) 


= Q(p, 01,02). Let 


, ^ A[Ei,E2,l3] ,, A;'|a;i,a;2l , ,, 

and B = , ,0. , A = —^^ and the map v 


{Yf - Yi) 


- xl) 


V : 


B- 


A' 


^ aiX2 
0103X2 


30 


E3F 


50 


E a/-| X-i 

+0212 2^ ^ 
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For this we do the following: 


LIB "GND.lib"; 

ring All = 0,(al,a2,a3,xl,x2,Y1,Y2,Y3),dp; 

int nra = 3; 

int nrx = 2; 

int nry = 3; 

ideal xid = xl~3-x2''2; 

ideal yid = Yl~3-Y2''3; 

ideal aid = a3~2+a3+l; 

poly y; 

int i; 

for(i=0;i<=30;i++) 

{ 

y = y + al^xl'i/factoriaKi); 

} 


//load the library 
//define the ring 
//number of a's 
//number of x's 
//number of Y's 
//define the ideal from A 
//define the ideal from B 
//define the ideal from k’ 


for(i=31;i<=50;i++) 

{ 

y = y + a2*x2*xl~i/factorial(i); 

} 

ideal f = al*x2,al*a3*x2,y; //define the map v 

desingularization(All, nra,nrx,nry,xid,yid,aid,f,"debug"); 


Example 24. We continue on the idea of Examples 141 [T6l The bordered matrix H 
defined above is equal to 


/ 2Fi + F2 W + 2^2 


H = 


V 


0 

1 


0 

0 

2 Y 4 


and hence G = N ■ G' is equal to 


G = Vl 


/ 0 0 

W + 2F2 0 

0 ¥4^ + 4V1V2 + 4 ^ 2 ^ 

y -2F4 0 


0 0 \ 

1 0 

0 0 

0 Fi + 21^ y 


yy + 4V4V2 + 4F2' 

-2F2 _ 51 ^ 1-2 _ 2 Y^ 
0 

3F1F4 


and s = 1. Using the definition of h in Equation [H we get that 


hi = 

h2 = 

h, = 


Fi - (0:2) • T3 - (01X2), 
F2- 


^3 


-|- (1\ 


a^xl + 2xl 

• Ti H---—— ■ T3 — (0103X2), 

203 + i 


0 

0 

0 




F1 + 2F2 y 


^3 - ( 3 ^ 2 ) • T2 - + ^“1^1 + ^“ 1^1 + + ]^aixl+ 

OiXi + Oi) , 


h/i — F4 + 


2x2 


(2O1O3 + Oi) 


■Ti- 


3 x 3 


O4 (203 + I)" 
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-n- 


Xn 


2 cIi (23 “h Qi 


•T4 + 


2 fli(Z 3 -|- Qi 












From Equation [2] we get that 


Qi — 
Q2 = 


XT, 


( 2 aia 3 + fli) 

4 X 2 


■T^ - 


3 x 1 


3alx\ + 3a3X^ + 3x^ 2 


cii ( 2 a 3 + 1 ) 

m2 12 X 0 

■ 2 


■ T1T3 + 


(2aici3 + cii) af (203 + 1) (203 + 1) 


2x^ 


■ T1T4 - 


3Xn 


.4 ■ ^ 1^3 - 
■T 3 T 4 


(2a3 + 1) 

9 x 2 


•T, 


3 5 


■tU 


(2ai(23 + cii) Ui (203 + 1) 

and therefore following the dehnition of g in Equation [3] we have 

gi = Qi + Ti + + dids + 04 ), 

92 = Q2 + T2. 

We print now the algorithm’s debug output using the line codes from Example 


This is the bordered matrix H: 
2*Y1+Y2,Y1+2*Y2,0,0, 


0, 

0, 

1,0, 

1, 

0, 

0,0, 

Z, 

2*Z, 

0,Y1+2*Y2 

This 

is G: 


0, 


0, 


G[l,3] 

Y1*Y4''2+2*Y2*Y4''2,0, G[2,3] 

0, G[3,2],0, 


0 , 

0 , 

0 , 


-2*Y4~3, 


0 , 


3=i=Y1*Y4''3,Y1*Y4''2+2*Y2*Y4~2 


G [1,3] =Y1~2*Y4''2+4*Y1*Y2*Y4''2+4*Y2''2*Y4~2 
G[2,3]=-2*Y1''2*Y4~2-5*Y1*Y2*Y4''2-2*Y2~2*Y4''2 
G [3,2] =Y1~2*Y4''2+4*Y1*Y2*Y4''2+4*Y2''2*Y4~2 


s = 1 
h = 

_[l]=Yl+(-x2~4)*T3+(-al*x2) 

_ [2] =Y2+ (-x2~3) / (2*al*a3+al) *T1+ (a3*x2''4+2*x2~4) / (2*a3+l) *T3+ 
(-al*a3*x2) 

_[3] =Y3+(-x2~4)*T2+(-al*xl''6-6*al*xl~5-30*al*xl''4-120*al*xl''3 
-360*al*xl''2-720*al*xl-720*al)/720 

_[4] =Y4+(2*x2~2)/(8*al''3*a3''3+12*al~3*a3~2+6*al''3*a3+al''3)*Tl+ 
(-3*x2~3)/ (8*al''2*a3~3+12*al~2*a3~2+6*al~2*a3+al~2) *T3+ 
(-x2~3)/(2*al*a3+al)*T4-1/(2*al*a3+al) 


m = 2 
QT = 

_[l] = (x2''2)/(4*al~2*a3''2+4*al''2*a3+al~2)*Tl''2+ 
(-3*x2~3)/(4*al*a3~2+4*al*a3+al)*Tl*T3+ 
(3*a3~2*x2~4+3*a3*x2''4+3*x2~4)/(4*a3''2+4*a3+l)*T3~2 
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_[2] = (-4*x2)/(16*al~4*a3''4+32*al''4*a3''3+24*al~4*a3~2+8*al''4*a3+al''4) 
*Tl~2+(12*x2''2)/ (16*al"3*a3"4+32*al''3*a3~3+24*al~3*a3~2+8*al''3*a3 
+al~3)*Tl*T3+(-9*x2~3)/(16*al~2*a3~4+32*al~2*a3''3+24*al"2*a3''2 
+8*al~2*a3+al~2)*T3~2+(2*x2~2)/(4*al''2*a3''2+4*al~2*a3+al~2) *T1*T4 
+(-3*x2~3)/(4*al*a3~2+4*al*a3+al)*T3*T4 

f = 

f[1]=Y1~2+Y1*Y2+Y2~2 
f[2]=Yl*Y4+2*Y2*Y4+(-x2~2) 

g = 

_[l] = (x2''2)/(4*al''2*a3''2+4*al~2*a3+al''2)*Tl''2+(-3*x2''3)/(4*al*a3''2+ 
4*al*a3+al) *Tl*T3+(3*a3~2*x2''4+3*a3*x2''4+3*x2~4) / (4*a3"2+4*a3+l) *T3~2 
+Tl+(al~2*a3''2+al~2*a3+al''2) 

_[2] = (-4*x2)/(16*al~4*a3''4+32*al''4*a3''3+24*al~4*a3~2+8*al''4*a3+al''4) 
*Tl~2+(12*x2''2)/ (16*al''3*a3''4+32*al''3*a3~3+24*al~3*a3~2+8*al''3* 
a3+al~3)*Tl*T3+(-9*x2~3)/(16*al~2*a3''4+32*al"2*a3"3+24*al''2*a3~2 
+8*al~2*a3+al~2)*T3~2+(2*x2''2)/(4*al''2*a3''2+4*al~2*a3+al~2) *T1*T4 
+(-3*x2~3)/(4*al*a3~2+4*al*a3+al)*T3*T4+T2 


Thus the General Neron Desingularization is a localization of D\Y,T]/{h,g) = 

D\T]Ka)- 


Example 25. In the case of Example 151 and ITTI we obtain that the bordered matrix 


/ 

3 ^ 1 ^ 

CO 

1 

0 

0 

\ 


0 

0 

1 

0 



1 

0 

0 

0 


V 

0 

-QY2Y4 

0 

-3E2" 

/ 


and hence G = N ■ G' is equal to 


0 \ 

0 

0 

31^2^ / 


/ 0 0 
-Wi 0 
0 Wi 0 

\ 6Y2Y4 0 - 18 Y^^Y 2 Y 4 - 
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G = Y^- 






and s = 1 . Using the definition of h in Eqnationdl we get that 
hi = Yi- (x®) • Ts + (-01X2), 

^2 = ^^2 + 2^ 2 ■ -1 ■ ^3 ~ (fllfl 32 ; 2 ) , 

Safog a. 


h. = Y. 




•To 




T^aiXi + ^aix^ + + \:aix\ + ^aixl + ^aixl + OiXi + Oi 

7 ! 0! 5 ! 4 ! 3 ! / 


/14 — ^4 + 


2xi 

,5^5 


■Ti 


2 x1 

3^5 


QagOg SafOg 

From Eqnation [ 2 ] we get that 


T 3 + 


^2 


2^2 


3 a|a 


T4 + 


2^2 • 


3 afa 


Qi — 


Q2 — 


^10 

X2 

27 afa^ 


■T^ 


X. 


12 


SOgh 


4^6 
3 


T^Tg + 


X, 


14 


2„6 


ata 


^ 6™16 

TlTg^ + “3^2 


X, 


16 


1 *^3 

,33 


OjO 


^2 to 2 2x^ 3aiaix^ - 3aiX^ 2 

■ Jg + -- — ■ 4 1-13 3-^- ■ -1^ 


3 afag 

20X3 

27 a|a^ 


.2-4 _ 


aittn 

2x11 


3 ai 75 ao 


2x1^ 

■ ^3 + ^ 

OgOo 


■ TgTg^ - 


3 ’ 

9 . 7.15 

2 ^^2 . j. 3 ^ 


Xo 


Oafog 

2x^ 

OgOg 


■ Tg 2 T 4 


•TgTg 


ig . ^U,g ^ 

9.7,12 .7,14 

4‘'4'2 'T' 'T' 'T' I "^2 


3^9 


oia 


3 afag 

33,8 

’ 4 X 2 _ ^2 

a?al ^ 


TgTgT4 + 

2a;^ 


afog 


■ 


Sa^al 


T1T4 + 


2x^ 

OgOg 


13 , 

4 

X2 

Safal 

•T3T4 


Tg^ + 


and therefore following the definition of 5^ in Eqnation [ 3 ] we have 

91 = Qi + Tg, 

92 = Q2 + T2 


To obtain this with Singular, we nse the same code lines as in Example |23l bnt 
we change the last one with 

desingularization(All, nra,nrx,nry,xid,yid,aid,f,"injective","debug"); 
Doing this, the algorithm will not compnte the kernel becanse of the injective 
argnment. 


Example 26. We do now the same compntations for Examples 151 [T 51 The bordered 


matrix H defined above is eqnal to 





( X2 

—Xg 0 

0 \ 


H = 

0 

0 1 

0 


1 

0 0 

0 



^ 0 

0 0 

—XgXg y 


and hence G = N ■ G' is eqnal to 




/ 

0 

0 

2 4 

X1X3 

0 

G = Yl- 

-Xgx| 

0 

0 

2 4 
X1X3 

XgX2x| 

0 

0 

0 






























and s = 1. Using the definition of h in Eqnationdl we get that 

hi = Yi + • T3 - (03X1), 

ha = la - {x\xl) ■ Ti + (x^XaX^) ■ T 3 - (a3a;a), 

h3 = U3 + • Ta - 

1 4 1 3 1 2 ^ 

+— 01 X 3 + ^012^3 + ^«i2^3 + «i 3^3 + ai I , 

h4 = I4 + • T4 + 1. 

From Eqnation [2] we get that 

Qi = 0 , 

Q2 = 0 

and therefore following the definition of g in Eqnation [3] we have 

9 1 = Ti 

92 = Ua. 

To compnte this with the library we do the following: 
ring All = 0,(al,a2,a3,xl,x2,x3,Yl,Y2,Y3),dp; 
int nra = 3; 
int nrx = 3; 
int nry = 3; 

ideal xid = x2~3-x3''2,xl~3-x3''2; 
ideal yid = Yl~3-Y2''3; 
ideal aid = a3~2-al*a2; 
poly y; 
int i; 

for(i=0;i<=30;i++) 

{ 

y = y + al*x3''i/factorial(i); 

} 

for(i=31;i<=50;i++) 

{ 

y = y + a2*x3''i/factorial(i); 

} 

ideal f = a3*xl,a3*x2,y; 

desingularization(All, nra,nrx,nry,xid,yid,aid,f,"debug"); 
The algorithm’s ontpnt is as expected: 

This is the nice bordered matrix H: 

(x2),(-xl),0,0, 

0 , 0, 1,0, 

1 , 0 , 0 , 0 , 

0, 0, 0,(-xl*x3~2) 

This is G: 

0, 0, (xl~2*x3''4)*Y4~2, 0, 
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(-xl*x3~4)*Y4~2,0, (xl*x2*x3~4)*Y4~2,0, 

0, (x1~2*x3''4)*Y4"2,0, 0, 

0, 0, 0, (-x1*x3''2)*Y4''2 


s = 1 
h = 

h[l]=Yl+(xl~3*x3''6)*T3+(-a3*xl) 

h[2]=Y2+(-xl''2*x3''6)*Tl+(xl~2*x2*x3''6)*T3+(-a3*x2) 
h[3] =Y3+(xl~3*x3''6) *T2+(-al*x3''7-7*al*x3~6-42*al*x3~5-210*al*x3''4 
-840*al*x3~3-2520*al*x3~2-5040*al*x3-5040*al)75040 
h [4] =Y4+ (-xl''2*x3''4) *T4+1 

m = 1 

QT = 

QT[1]=0 
QT[2]=0 
f = 

f[l]=(x2)*Yl+(-xl)*Y2 
f [2] = (xl*x3~2)*Y4+(-xl*x3''2) 


g 

_[1]=T1 
_[2]=T2 

Thus the General Neron Desingularization is a localization of D\Y,T^,T/^/{h) = 

D[n,T,]. 


Example 27. We do now the same computations for Example fT^ In this example, 
the computations are much more complicated. The output is unfortunately too big 
but we will try to describe the result. 

The bordered matrix H dehned above is equal to 




/ 

0 X 2 -Y 3 

xs ■ F2 -2 ■ n 

0 \ 





xj-Y 2 xl-Yi 

- 2-^3 0 

0 



E = 


0 1 

0 0 

0 





1 0 

0 0 

0 




V 

0 0 

4Y4Y, 4Y3Y, 

4Y3Y, J 


and hence G = N ■ G' 

is 

equal to 





( 

0 


0 

0 

16 Y 3 %^ 



0 


0 

16 Y 3 %^ 

0 


G = Yi- 


0 



8 xf ■ FiFsn' 

8 x! ■ Y 2 Y 3 Yi 



-%YiY 4 


- 4 a ;2 • EaWn 

G[ 4 , 3 ] 

4 xfx 2 ■ YiY 3 Yi 


V 

8E3W5 0:2 

Y^Y^^Y, + 2^42^5 

G[ 5 , 3 ] 

G[ 5 , 4 ] 
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0 \ 
0 
0 
0 


4V3V4 / 






where 


^[4, 3] = 4x^,X2 • + 2x2 • n, 

^[ 5 , 3 ] = - 4 x?X 2 • Y,Y2Y3Y5 - 2 x 2 • - 2 x 2 • Y^Y^^ and 

^[ 5 , 4 ] = - 4 x?X 2 ■ YiYsY^ - 2 x? ■ Y2Y/Y3 


and 


s = 0*02X22 — 20^03X2* + 0202X2"^ — 20^03X22 + 20^02X2* — 0^02X22 + 20103X2*+ 
20*02X3° — 4 ofo|x 2 ^ + 40203X32 — 4 o° 02 X 2 ° + 40^02X2* + 03X32 + 20103X2^+ 
0*X2 — 20^03X3 + 30^03X3° — 2o°02X* + 20503X3 + 05O2X* — 20^02X2 + 
20io5x2 + 20503X3 + 2o5x2 — 20^X3 + 20io5x2 + 20^02X3 + 1 


Using the dehnition of h in Equation [H we get that 


/h*\ 


( Yi-y'i \ 


/ Ti \ 

h 2 


Y 2 -y '2 


T2 

hs 

= s ■ 

y^-y's 

- xlG{y') ■ 

n 

/14 


Y,-y', 


T4 

\h ) 


\y^-y', 


\T 5 J 


where 


y[ 


y'2 


y's 

y'i 

2/5 


— 180503X32 + 50503X3* + 703X32 — 180103X3* + 15 o 5 o|x 2 ° — 4 o 5 o 5 x 2 — 
603X3° + 150103X3 — 12o5o5x5 + 3o5o5x 2 + 503x5 — 12oio5x5 + 9o5o5x5— 
2 o 5 o 2 x 5 — 4o5x5 + 9oio5x5 — 6 o 5 o 2 x| + o5x5 + 3o5x5 — 601O2X5 + 3o5x5 — 

202X3 + 30 iX 2 + 1 


a}^x^ + 2o5°02x 52 — o 5 *x 5 * + 050^x52 — 2 o 5 o 2 x 5 * + o 5 °x 5 ° — o 5 o 5 x 5 * + 
2o5o2x 5° — 0^x5 + o5o5x5° — 2 o5o2x 5 + o5x5 — ofo^X^ + 2 o5o2x 5 — o5x5+ 


o^o^x^ — 2 a\a 2 x\ + o 5 x 5 — o 5 o 5 x 5 + 20 ^ 02 x 5 — ofx^ + o 5 o 5 x 5 — 2 o 5 o 2 x 5 + 

0 ^x 5 — OiO^X^ + 2 o5o2x| — o 5 x 5 + O^x^ — 20 i02x5 + o 5 x 5 + 2 o2x5 — O1X2 + 


1 


Oix 5 + X2 
02 x 5 + X2 

o 5 4 o 5 

~r ^2 - 

4 2 


+ 0102 
"1 


x 5 + 


^2 9 

-^- ^2 



However, the output is too big to be printed. Following the idea in the above ex¬ 
amples, we compute Q and g. This is even bigger than h so we print the numerators 
and denominators of the coefficients just till degree 10 in the xds. However in some 
cases even the terms till degree 10 will be too many to write down and hence we 
will print just the hrst terms and . 

As a small remark, Q3 contains also terms in degree 3 in the Tj but the numerator 
of the coefficients have power greater than 10 and therefore they do not appear in 
our shortcutting. 
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Qi 


Q2 


Qs 


3 a‘fx‘lx^— 2 aix^X 2 -\- 4 :a 2 X^x^-\-x^x^ 


•T1T4 


4 aia 2 ^|+ 8 aia 2 a;^+ 4 aia; 2 H- 4 a 2 X 2 + 8 a 2 a: 2+4 4 a 2 a; 2 + 8 a 2 X 2+4 

^ —afx 2 ^-\-a^x^— 2 a^a 2 X 2 ^—a^x^-\- 2 aia 2 X^^aiX 2 — 0 ' 2 ^ 2 ^—^ 0 ’ 2 X 2~^2 

4 ai a 2 ^ 2 +^^i ^ 2 ^ 24 "^^! ^ 2 + 4 a 2 X 2 + 8 a 2 0:2+4 
—5a|x2*^H-4aJa:| —12aja2o:2*^—3a^a:|+80.10.20:2+2ciia:^—6a2a:2'^—4a2o:2—o:| 
16 a^ a 2 o:|+ 32 a^ 02 o:^+ 16 a^ a:^+ 32 ai 020:2+640102 o:|+ 32 oi 0:2+16020:2+3202 0:2+16 

^ o^o:Jx|+ 2 oJx 2 *^+ 2 oio:jO: 2 + 4 oix|+o:jO:|+ 2 x| 
401020:2+801020:2+4010:2+4020:2+8020:2+4 
^ 010:^0:2—2010:2+2o2o:^x|— 4020:2^+o:Jo:|— 2 o:| 

80^0^0:2 + 160^ 020:2+80^X2 + 1601020:2+3201020:2+16010:2+8020:2+16020:2+8 
— x|x2+ 4 x^x 1 — 4 X 2 *^ 

16 o^ 03X2+320^ 02X2+160^X2+320102X2+640102X2+3201x2+1602X2+3202X2+16 

^ 6 oJxJx 2 — 3 oix^x^ + 6 o 2 xJx|+x^x| 

80^02x2 + 160^ 02x2+80^x2 + 160102x2+320102x2+1601x2+802x2+1602x2+8 

_x^ “ 1 “ 2x^ x^ 

80^02X2 + 160^ 02X2+80^X2 + 160102X2+320102X2+1601X2+802X2+1602X2+8 


+ 


16 o^ 03X2+320^ 02X2+160^X2+320102X2+640102X2+3201x2+1602X2+3202X2+16 


40^X2+801X2+4 

O^x^xl—OlX^X 2 + 202 X^X 2 +X^X 2 


3 a 2 a; 2 a; 8 + 3 aix 2 x 7 - 2 a 2 a; 2 x 8 +x 2 x 6 
• H-—- ' 2-^3 + 


20^X2+401X2+2 


20^X2+401X2+2 

rri rp , X^X 2 -|-^X;l^X 2 

-I2-I4 + 


40^X2+801X2+4 


T3T4 + 


+ 


+ 


+ 


20^X2+4 oiX 2+2 40^X2+801x2+4 

20^X2—40^02X2*^—20^X2+20102X2+201X2—202X2—2x2 
40102X2+I20102X2+I20102X2+401X2+402X2+I202X2+1202X2+4 
Tojx®-28oJ 02X2*^—70^X2+210102X2+701X2—2102X2*^—2102x1—7x2 
...+48ojX2+48oi 03X2+1440102X2+I440102X2+4801X2+I602X2+4802X2+4802X2+I6 
^ 2oJx^X2+2o^X2*’+4oiX^X2+4oiX2—202 xJx 2 —202X2*^ + 2 XjX 2 + 2X2 
4oi 03X3+12oi 03X2+12oi 02X3+401X2+402X2+I202X2+I202X2+4 
70 jX^X 2+4 o^X 2'^ + 14 oiXjX 2+8 oiX 2+702 XjX 2+402X2'^ + 7 x^X 2+4X2 
...+240^x2+240102x2+720102x2+720102x2+2401x2+802x2+2402x2+2402x2+8 
— 7 x|x 2—8x^X3—4X2*^ 

...+48ojX2+48oi 03X2+1440103X2+1440102X3+4801X2+I602X2+4802X2+4802X2+I6 

6o^XjX2—4aiX^X2+602XjX2+2x^X3 
4oi 03X2+120103X2+120102X3+401x2+403X3+1203X2+1202X3+4 
21 o^xJx2 — 1401 X^X 3 + 3502 X^X 3 + 7 x^X 3 

...+24 oJx 2+24 oio|x^+ 72 oi 02 x|+ 72 oi 02x^+2401x2+80^x^+2403X2+2402X2+8 

_7X^X3_4X^X3 

...+ 24 a^x|+ 24 aia 2 X 2 + 72 aia|x 2 + 72 aia 2 a: 2 + 24 aiX 2 + 8 a 2 X 2 + 24 a|x 2 + 24 a 2 X 2+8 

20 


+ 


Tf 

T1T2 

-^2 

T1T3 

T2T3 

rp2 

-‘-3 

T2T4 

T3T4 

T1T2 

rri2 

-^2 

T1T3 

T2T3 

-^3 

T1T4 

T2T4 

T3T4 




The General Neron Desingularization is a localization of D\Y,T]/{h, g). For this 
example we will need a function 

invpCpoly p, int bound,string param,string variab) 

which computes computes the inverse of p till order bound in Q(parain)[variab]. 
The input for this example is the following: 

ring All = 0,(al,a2,xl,x2,Yl,Y2,Y3,Y4),dp; 

int nra = 2; 

int nrx = 2; 

int nry = 4; 

ideal xid = xl''2-x2~3; 

ideal yid = Y3''2-xl''2*Yl*Y2,Y4''2-x2*Y2*Y3; 

ideal aid = 0; 

poly yl,y2,y3,y4; 

y3 = l+al*x2; 

y4 = l+a2*x2''2; 

string as,xs; 

if(nra != 0) 

{ 

as = stringCvar(1)); 
for( int i=2;i<=nra;i++) 

{ 

as = as+","+string(var(i)); 

} 

} 
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if(nrx!=0) 

{ 

xs = stringCvar(nra+1)); 
for(int i=nra+2;i<=nra+nrx;i++) 

xs = XS+" ,"+string(var(i)); 

} 

} 

yl = y3~3*invp(y4''2,12,as,xs); 
y2 = y4~2*invp(y3,12,as,xs); 
y3 = x2*y3; 
y4 = x2*y4; 

ideal f = yl,y2,y3,y4; 

desingularization(All, nra,nrx,nry,xid,yid,aid,f,"injective","debug"); 

Remark 28. Our algorithm works mainly for local domains of dimension one. If 
A' is not a domain but a Cohen-Macaulay ring of dimension one then we can build 
an algorithm in the idea of the proof of Theorem [201 In this case it is necessary 
to change B by an Elkik’s trick | 3 ] (see in Lemma 3 . 4 ], [131 Proposition 4 . 6 ], [TTl 
Corollary 5 . 10 ]). The algorithm and as well Theorem 1201 might be also build when 
A' is not Cohen-Macaulay substituting in the proofs d by a certain power (B such 
that (0 :a d^) = (0 :a Such algorithm could be too complicated to work really. 

On the other hand, if we restrict our present algorithm to the case when A' is the 
completion of A then we might get a faster algorithm using the idea of the proof of 
Theorem [201 This algorithm could be useful in the arc frame. 
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